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We present a functional renormalization group calculation of the effect of strong interactions on 
the shape of the Fermi surface of weakly coupled metallic chains. In the regime where the bare 
interchain hopping is small, we show that scattering processes involving large momentum transfers 
perpendicular to the chains can completely destroy the warping of the true Fermi surface, leading 
to a confined state where the renormalized interchain hopping vanishes and a coherent motion 
perpendicular to the chains is impossible. 
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Introduction. Electron-electron interactions can 
strongly modify the Fermi surface (FS) of a metal. A well 
known example is the Pomeranchuk transition, where 
the symmetry of the FS is spontaneously broken due to 
strong interactions in the zero-sound channel How- 
ever, there are other quantum phase transitions associ- 
ated with the geometry or the topology of the FS with- 
out symmetry breaking, such as the Lifshitz transition 
0, S] or the truncation transition 0, 0], where certain 
sectors of the FS are washed out by interactions, while 
others remain intact. Another example is the interaction- 
induced confinement transition, which has been proposed 
by Clarke, Strong, and Anderson more than ten years 
ago Q : they considered metallic chains with small inter- 
chain hopping t±,o. For weak interactions, the FS con- 
sists then of two disconnected weakly curved sheets as 
shown Fig. [T] The amplitude of the warping of the FS is 
proportional to the renormalized interchain hopping t±. 
Clarke et al. Q suggested that, at least for sufficiently 
strong interactions and small i±,0j the renormalized t± 
vanishes. The true FS is then completely flat, so that a 
coherent motion of the electrons in the direction perpen- 
dicular to the chains is not possible (confined coherence) . 

In the past decade the confinement problem in weakly 
coupled metallic chains has been studied by many au- 
thors @, 0, H, @, EH , but the results have not converged 
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FIG. 1: (Color online) FS of a two-dimensional array of 
weakly coupled metallic chains. The dashed lines mark the 
boundary of the first Brillouin zone in the direction perpen- 
dicular to the chains. 



due to a lack of controlled methods. A simple one-loop 
calculation 0, Q suggests that the renormalized inter- 
chain hopping vanishes if the anomalous dimension r\ of 
the Luttinger liquid state without interchain hopping is 
larger than unity. However, this argument does not take 
the renormalization of r\ by interchain hopping into ac- 
count. Indeed, more refined calculations by Arrigoni (loj 
suggest that higher order corrections are important and 
possibly lead to a finite t± even for rj > 1. In this 
Letter we shall re-examine this problem using a novel 
functional renormalization group (RG) approach involv- 
ing both fermionic and bosonic fields [ll|, [l2| • Our main 
result is that the regime of confined coherence proposed 
in Ref. [6] indeed exists, so that strong interactions can 
give rise to a non-Fermi liquid normal state in quasi-one- 
dimensional metals. 

Model. We start from an effective low energy model for 
spinless fermions with linearized energy dispersion and 
density-density interactions. The Euclidean action is 

S[^,ip] = ^2 [-iv + av F 5k\i + fJ.o(k±)]i>K a 'ipK a 
a Jk 

+ f aa 'PRaPK a > , (!) 

aa' ^ 

where fJ.o{k±_) = — T>(kp,i0) is a counter-term involv- 
ing the exact self-energy at the true FS k — kp and 
zero frequency, and 5k\\ = fcii — akp(k±). Here fcii is 
the component of the two-dimensional lattice momen- 
tum k parallel to the chain direction, and kj_ is the cor- 
responding perpendicular component. The FS kp can 
then be parameterized as fcy = akp(kj_), where a = ±1 
labels the two disconnected sheets of the FS, see Fig. [T] 
We neglect the /c^-dependence of the Fermi velocity v F . 
The chiral fields ipKa arc defined in terms of the usual 
Fermi fields ipkpkx,^ via V>K« = ip a k F (k ± )+Sk h k ± ^, and 
the chiral densities p Ra are p Ra = J K tpKa^K+Ra- We 
use collective labels K — (5ku, k±,iu) for fermionic and 
K = (k^kj_,iu)) for bosonic fields, where to and u are 
Matsubara frequencies. The integration symbols are 
f f f A ll dSk \\ f 1 r _ r r^ii d>H r da 

JK Jk ± J-Ai| 2tt J 2tt allu JK Jk ± J —Am 2tt J 2tt 
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where for later convenience we have introduced the no- 
tation j k± = /_ A ^ and f k± = §L. Here A„ 

is a bandwidth cutoff, Aj_ = tt /a± is the width of the 
Brillouin zone in transverse direction, and A^ and A|| re- 
strict the momentum transfered by the interaction in the 
directions parallel and perpendicular to the chains. We 
assume that A|| <C min{fcF (&_!_)}, so that the interaction 
f aa i in Eq. (fTJ) does not transfer momentum between the 
two disconnected sheets of the FS. However, the trans- 
verse momentum transfer cutoff Aj_ can be of the order 
of the transverse width A±_ — n/aj_ oi the Brillouin zone, 
so that transverse Umklapp scattering is possible. For 
simplicity we set Aj_ = A± and call A|| — Aq. 

Self- consistent perturbation theory. To begin with, let 
us calculate the FS within second order self-consistent 
perturbation theory. Using the procedure outlined in 
Ref. 13|, we obtain the following integral equation for 
the difference 6kp(k±) = kp(k±) — kpfl(k±) between 
the true Fermi momentum kp(k±) and the corresponding 
kF,o(k±) without interactions at the same density, 



5kp(k±) 
A 
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J(A(fc ± ,fc ± );A(fcl,fc ± )) , (2) 



where A(fcj_, k±) = [kp(k±) — kp(k± + k±)]/Ao, and 



J(A;A') 



A' 



In 



'\2 



4 _ (A - A') 
(A + A') 2 



(3) 



The dimensionless couplings 52 and 34 are defined via 
271-54 = Vof++ = v Q f — and 2irg 2 = vof+- = fof-+, 
where the factor i/q = A^^vp)^ 1 = (a^vp)^ 1 is in- 
troduced for convenience. We have solved Eq. ([2]) nu- 
merically, but for small ij_ we can also obtain an ap- 
proximate analytic solution using the fact that in this 
case the dominant renormalization of the FS is due to 
the logarithmic term in Eq. ([3]). Suppose that the bare 
FS is of the form kFfi(k±) = kp + to cos(k±a±) where 
to = 2t±^/vF <C Ao and the average kp is fixed by 
the total density. The renormalized FS is then given by 
kp(kj_) = kp + 1 cos(fcj_aj_) + . . ., where t is proportional 
to the renormalized nearest neighbor interchain hopping, 
and the ellipsis denotes higher harmonics corresponding 
to longer range hoppings. From the numerical solution 
of the integral equation @ we find that for 32,54 -C 1 
the higher harmonics are indeed small. Then Eq. @ 
can be reduced to a transcendental equation for t, which 
to leading logarithmic order in t/Ao can be written as 
t/t = [l+-R(t)] -1 s with 



(4) 



A similar relation has been obtained previously [HI, [3| 
for the difference between the Fermi momenta associ- 
ated with the bonding and the anti-bonding band in two 



coupled spinless chains. Note that to first order in the 
bare interaction R(t) oc 174, so that a repulsive interaction 
#4 > reduces the warping of the FS, while for g± < 
the warping of the FS is enhanced. However, for t±_ — > 
the logarithmic term proportional to g\ always dominates 
and predicts an interaction-induced reduction of the FS 
warping, irrespective of the sign of the interaction. 

Functional RG approach. We now generalize the RG 
approach developed in Ref. [13] in the context of a sim- 
plified two-chain model to the more interesting two- 
dimensional case considered here. The method has been 
described in detail previously [l2j], so that we will be 
rather brief here. In the momentum transfer cutoff 



scheme H( we decouple the density-density interaction 
by means of a bosonic Hubbard-Stratonovich transfor- 
mation and then use the maximal momentum carried by 
the boson field as flow parameter A of the RG. Our ini- 
tial cutoff is thus A = Ao = Am. Eliminating boson 
fields with momenta in the range A < |fcii| < Ao we 
obtain a new effective action, whose vertices are deter- 
mined by a formally exact hierarchy of functional RG 
flow equations. To calculate the true FS, we need the 
flow equation for the relevant part rj(fcx) of the irre- 
ducible fermionic self-energy H\(K,a), which is defined 
via n(fcx) = Zi(k±)[Y,A(kp, iO, a) + jj,o(k±)]/vpA. Here 
I = ln(Ao/A) and Zi(k±) is the flowing wave- function 
renormalization. The functional RG flow equation for 
fi(k±) is of the form 



din(k±) = [1 - vi(k±)}n(k x ) + ti(k±) , 



(•5) 



where T)i(k±_) — —dihxZi(Jc±) is the flowing anomalous 
dimension. An approximate expression for the inhomo- 
geneity r^fcjj is given in Eq. {2} below. As long as 
Voc(k±) < I, the FS is well defined. The shift 5kp{k±_) 
of the FS due to interactions can then be obtained from 
the requirement that the initial value ro(fcjJ should be 
fine tuned so that the relevant coupling ri(k±) flows into 
a fixed point [3] . This leads to the following exact inte- 
gral equation for the FS, 



Skp(k±) 
A 



ro(kA 



dle- l+ ^ dtr >^ kx) ti(k ± ) 



Using the same truncation strategy as in Ref. [121] . we 
approximate 



dqde S(\q\ - l)[F l (q,ie,k±)] aa e 



(TO 



X7,(fcj 



(2tt)2 
k±)ji{k A 



aq - Ai(k ± ,k A 
a k 1 ) , 



(7) 



where A;(A;j_,fcx) = kpj(k±) — kp^i{k±_ + fcj_) with 
kp,i(kx) — kp{kx)/A - n(kj_). The inverse of the 2x2- 
matrix Fi(q, il, k±) is defined via 

\Fl{q,ie,k±)]~a> = [vof\Ztf +S aa >fli{q,ie,k±,a) , (8) 
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FIG. 2: Diagrammatic representation of the flow equatii 
(jlll) for the three-legged vertex with one bosonic (wavy lin 
and two fermionic (solid lines with arrows) external lef 
The thick wavy line with a slash denotes the bosonic sinj 
scale propagator. Additional contributions involving irrel 
vant higher order vertices are omitted, see Refs. [TTI. fl^] ■ 



where f is a matrix in chirality space with elcmcn 
[f]oa' — faa', and IT/ (q, ie, k± , a) is the rescaled pole 
ization associated with fermions of chirality a, for whii 
we use the adiabatic approximation [r3 | 



TIi(q,ic,k±,a) = — 



A t (k±, k±) + aq 



2tt J k± Ai(k±,k±) + aq-ie 
x ji(k±, k±)ji(k± + fcj_, -kj_) • ( 

The anomalous dimension is in this approximation 

dqde S(\q\ - l)[F l (q,ie,k±)] aa 



k± j (2tt) 2 [ie-aq-K^k^k^Y 
xyi(k±,k_ L )yi(k± + k±,-k±) . (1 

Finally, the dimcnsionlcss vertex "fi(k±,k±) with o: 
bosonic and two fermionic external legs (where k± labf 
the incoming fermion and k± labels the boson) satisfi 
the flow equation 



dm(k ± ,kj_) = --[r]i(k±) +m(k±+ k±)]'yi(k±,k±) 



dqde 



S(|«|-l)[F,(g,ie,fcl)] c 



li(k± + k' ± ,k±)'yi(k±, k' ± hi(k± + k ± + k' ± ,-k' ± ) 
[ie- aq ~ A t (k±, k' ± )] [ie - aq - A/(fcj_ + k±,k' ± )] 

(11) 



with initial condition "fi=o(kj_, k±) — 1. A graphical rep- 
resentation of Eq. (JTTJ) is shown in Fig. [5] 

Results. Eqs. ©-(dU form a closed system of flow 
equations for the rescaled self-energy at the FS ri(k±), 
the flowing anomalous dimension r]i(k±), and the three- 
legged vertex yi(k±, k±). Of course, these equations can 
only be solved numerically, but the qualitative behavior 
of the solutions can also be extracted analytically. To 
begin with, let us establish the relation with the per- 
turbative Eq. ([2]). We set g± — from now on, be- 
cause the dominant renormalization of the FS is due to 
the cfc-process. In the simplest approximation, we set 
ji(k±,k±) ks 1, rji(k±) w and replace the flowing FS 
k F ,i(k±) = k F (k^)/A - n(k±) by k F (k±)/A. Then we 
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FIG. 3: (Color online) Typical evolution of the vertex 
7;(fcx,fc±) for different values of the flow parameter /. To 
evaluate the flow we have expanded in Eqs. (J7J) and pip 
up to g%. We have assumed a harmonic bare FS with am- 
plitude to/Ao — 10 -3 and bare coupling gi = 0.4 From 
left to right / = \l c ,l c ,Zl c and 7 min = 0.999989,0.981,0.85, 
where the crossover scale l c (see text) can be approximated 
by l c m — ln(2£o/Ao) for small gi- 



obtain from Eqs. © and to leading logarithmic order 
= [(l-sIrM / Mk±,k±)hx\A(k ± ,k±)\ ■ 

A o Jk ± 

(12) 

Expanding in harmonics we obtain as before t/to = [1 + 
iZ(t)] -1 , but now with 

R(t) « [(l-.g 2 2 )^-l]lu(A /|i|), (13) 

which reduces to Eq. (01 to leading order in g 2 - Obvi- 
ously, R{t) — > oo for <j 2 — > 1, indicating a confinement 
transition at strong coupling, where the renormalized in- 
terchain hopping t vanishes. However, from our previous 
work on two coupled chains [i~2l ] we know that vertex cor- 
rections and wave-function renormalizations can possibly 
change this scenario. 

To investigate this, let us first consider the RG flow 
of the vertex 7z(/c_i_, k±) numerically. Representative re- 
sults are shown in Fig. [3] Obviously, the dependence of 
"fi(kj_, k±) on the fermionic momentum fcj_, which devel- 
ops at intermediate scales /, is smoothed out again at 
the scale l c where the reduced cutoff A c = A e~ lc be- 
comes comparable with the warping of the renormalized 
FS. Defining the flowing dimensionless nearest neighbor 
interchain hopping ti via k F ,i = k F /A+ti cos(fcj_a^) + . . ., 
the crossover scale l c can be defined self-consistently via 
2i;=; c = 1. Because of the weak dependence of 7/ (k±, k±) 
on the fermionic momentum fcj_, we may approximate 
-fi(k ± ,k ± ) w 7;(0,fcj_) = ji(k±). Moreover, using the 
fact that close to the confinement transition |t;| -c 1, 
we obtain from Eqs. (|5jl- (|lip to leading order in £/ the 
following RG flow equation for the effective interaction 



9i(k±) = 92li{k± 



9i9i{k±) 



4sm 2 (k ± a ± /2)g l {k ± )uft 2 



«f[l 



0" 



,213 



(14) 



with ui = gi (ir/a±). Note that the flow of gi [k±) is driven 
by the component gi(n/a±) of the interaction involving 
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FIG. 4: (Color online) Renormalized nearest neighbor in- 
terchain hopping t± and renormalized interaction u = 
<?oo(7r/ax) as a function of the bare interaction 32 for to = 
2£_l,o/(«fAo) = 0.1 as obtained from the numerical solution 
of Eqs. (0 and For the perturbative curves we have 

expanded in these equations up to order g\. 











de-confined phase 


/ confined phase 



FIG. 5: (Color online) Projection of the RG flow in the t- 
u plane. The trajectories are obtained from the numerical 
solution of Eq. (|14p and ti = (l — fji)ti, where f\i is the weighted 
FS average of the flowing anomalous dimension r)i(k±) defined 
in the text. 



momentum transfer fcj_ = ir/a±; in a simplified two-chain 
model this corresponds to the pair-tunneling process [3] • 
The flow of the rescaled interchain hopping ti is deter- 
mined by diU = (l — fji)ti + 0(tf), where fji is the weighted 
FS average fji = 2 J k ^ sm 2 (k±a±/2)r]i(k±) of the flowing 
anomalous dimension, which for small i; can be approxi- 
mated by rji(k±) = [1 — gf(k±)]~? — 1. From the numer- 
ical solution of these equations we find that, for a given 
t± t o <C VpAo, there exists a critical value of the bare 
interaction (72 where the renormalized t = vanishes, 
corresponding to a confinement transition to a state with 
vanishing interchain hopping t± and flat FS. In Fig. [4] we 
show the ratio tx/tj_,o together with the renormalized 
interaction u — g QO (Tr/a±) as a function of the bare in- 
teraction for small to = 2t±^/(vF A ). A projection of 
the RG flow in the t-u plane is shown in Fig. [5l At the 
confinement transition the system is certainly not a Fermi 
liquid, because the anomalous dimension r]i(k±)js larger 
than unity and a sharp FS cannot be defined [la ]. The 
physical properties of the model in the confined phase 
remain to be explored. 

Summary and conclusions. We have shown that suffi- 
ciently strong interactions involving momentum transfers 
parallel to the chains can lead to a confinement transition 
in highly anisotropic quasi-one-dimensional metals. At 
the confinement transition the curvature of the two dis- 



connected sheets of the FS vanishes, so that a coherent 
motion of the electrons perpendicular to the chains is not 
possible and the electronic motion is one-dimensional. 
Our calculation thus supports the existence of a confined 
state in quasi-one-dimensional metals, as originally sug- 
gested by Clarke, Strong, and Anderson [6]. In the con- 
fined state the system is a non-Fermi liquid with large 
anomalous dimension. From the numerical solution of 
our functional RG equations we found no evidence for a 
truncation transition 0, [1| , where only certain sectors of 
the FS are washed out by interactions. 

Because in this work we have considered only spin- 
less fermions and interactions which do not transfer mo- 
mentum between the two disconnected sheets of the FS, 
our model ((TJ) is too simple to describe the competi- 
tion between confinement and the tendency to develop 
some kind of long-range order, such as charge-density or 
spin-density waves. In principle, spin fluctuations can be 
taken into account with the help of another Hubbard- 
Stratonovich field, but the analysis of the resulting RG 
equations for the coupled boson-fermion model requires 
a substantial extension of our calculation. 

We thank A. Ferraz and F. H. L. Essler for useful dis- 
cussions. 
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